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ABSTRACT 

During the late stages of stellar evolution in massive stars (C fusion and later), the fusion 
luminosity in the core of the star exceeds the star's Eddington luminosity. This can drive vig- 
orous convective motions which in turn excite internal gravity waves. The local wave energy 
flux excited by convection is itself well above Eddington during the last few years in the life 
of the star. We suggest that an interesting fraction of the energy in gravity waves can, in some 
cases, convert into sound waves as the gravity waves propagate (tunnel) towards the stellar 
surface. The subsequent dissipation of the sound waves can unbind up to several Mq of the 
stellar envelope. This wave-driven mass loss can explain the existence of extremely large stel- 
lar mass loss rates just prior to core-collapse, which are inferred via circumstellar interaction 
in some core-collapse supemovae (e.g., SNe 2006gy and PTE 09uj, and even Type Iln su- 
pernovae more generally). An outstanding question is understanding what stellar parameters 
(mass, rotation, metallicity, age) are the most susceptible to wave-driven mass loss. This de- 
pends on the precise internal structure of massive stars and the power-spectrum of internal 
gravity waves excited by stellar convection. 
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1 INTRODUCTION 

Many massive stars appear to lose a significant fraction of their 
mass in episodic outbursts rather than continuous line-driven winds 
(e.g., |Bouret et al.|2005[ [Smith & Owocki|2006| . There is strong 
evidence from observations of luminous supemovae (SNe) that 
the most extreme version of this phenomena is the ejection of 
~ 1 - IOMq of the stellar envelope in the last year to decade of 
massive stellar evolution (in a very small fraction of massive stars). 
In particular, the interaction between an outgoing supernova shock 
and such ejecta can explain some of the most optically luminous 
SNe yet detected (e.g.,|Smith & M cCray|2007| l, including, e.g., SN 
2006gy ( [Smith et al.]2 007 , Ofek et al .|2007[ > and perhaps the emerg- 
ing class of hydrogen-poor ultraluminous SNe jQuimby et al.|20lT| 
[Chomiuk et al.[20TT) . In several cases, the late-time light curve dis- 
favors one alternative explanation, that the luminosity is powered 
by unusually large amounts of radioactive Ni and Co (e.g.,|Miller| 
[et al.|20 10] [Chomiuk et al.|201ip . Related evidence for prodigious 
mass loss in the last few years of stellar evolution comes from SNe 
like PTF 09uj, which was interpreted as shock-breakout from an 
extremely dense circumstellar wind (Ofek et al .|2010^ . Even more 
typical Type Iln SNe (i.e., those displaying narrow emission lines 
indicative of circumstellar interaction) appear to require progeni- 
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tor mass loss rates exceeding ~ lO^'Moyr"' Kiewe et al.|2012 i, 
far larger than can be explained by continuous mass loss processes 
operating in massive stars. 

One of the central puzzles in the circumstellar interaction sce- 
nario for these SNe is why a massive star should lose a substantial 
fraction of its mass in only ~ 10"*"^ of its lifetime. This is required 
in order for the supernova shock to encounter the stellar ejecta at 
radii ~ 100 AU where the shock is particularly radiatively efficient. 

The most important change that occurs in the late stages of 
massive stellar evolution is the onset of prodigious neutrino cool- 
ing in the core of the star, associated with the high temperatures 
required for carbon fusion and beyond (see, e.g., [Woosley et al.[ 
[2002[ for a review). The fusion luminosities exceed Eddington for 
C fusion and later and become particularly large during the last 
year in the life of the star. The rapid fusion and neutrino cool- 
ing in turn accelerate the nuclear burning so that the characteristic 
nuclear timescale t„„c (analogous to the main sequence timescale 
for H-buming) decreases to a ~ yr for Ne and O fusion and a ~ 
day for Si fusion (with the exact values depending on stellar mass, 
metallicity, and rotation). Because neutrino cooling and fusion have 
different temperature sensitivities, it is in general not possible for 
neutrino cooling to balance nuclear heating everywhere in the core 
of the star, although in a volume averaged sense the two balance. 
This local dilference between heating and cooling drives convec- 
tion which carries a significant convective luminosity ~ lOi % of 
the fusion luminosity. 
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Table 1. Late Stages of Massive Stellar Evolution 



Stage Duration (t„„c) Lfusio„ (Lg) Mach (Alconv) (s) 

Carbon ~ 10^ yr ~ 10** ~ 0.003 ~ lO'*"' 

Neon ~ 1 yr ~ lO' ~ 0.01 ~ 10^ 

Oxygen ~ 1 yr ~ lO'" ~ 0.02 ~ 10^ 

Silicon ~ 1 day ~ lO'^ ~ 0.05 ~ 10^ 



Fusion luminosities, durations, convective Mach numbers, and convective 
turnover times for core fusion during the late stages of stellar evolution of a 
~ 25Mq star (based on |Woosley et al.|2002^ . Precise values depend some- 
what on stellar mass, metallicity, and rotation. Convective Mach number is 
an order of magnitude estimate assuming that a significant fraction of the 
fusion luminosity is locally canied by convection; depending on the stel- 
lar parameters, core carbon fusion may not be convectively unstable. Shell 
fusion of C. O, Ne. etc. can produce similarly vigorous convection. 

In this Letter we argue that wave excitation by vigorous con- 
vection in the late stages of stellar evolution is capable of driving 
the strong mass loss suggested by circumstellar interaction in lumi- 
nous core-collapse SNe. Such wave excitation has been explicitly 
seen in numerical simulations of carbon and oxygen shell fusion by 
Amett and collaborators (e.g., |Meakin & Arnett|2006j ); there is also 
closely related numerical work in the context of solar convection 
(e.g., [Rogers & Glatzmaier|2005|[2006) . Here we provide analytic 
estimates of wave excitation in evolved massive stars and discuss 
the resulting implications for mass loss in the last ~ year of stellar 
evolution. We begin by summarizing some of the key properties of 
convection during carbon fusion and later and the excitation of in- 
ternal gravity waves by such convection (^. We then calculate the 
conditions under which a super-Eddington flux of waves excited in 
the core of the star can tunnel through to the stellar envelope (^. 
Finally, in ^we discuss the implications of our results and direc- 
tions for future work. 

Because there are significant uncertainties in the interior struc- 
ture of massive stars during the evolutionary phases of interest, we 
focus on elucidating the general conditions required for efficient 
wave-driven mass loss. We defer to future work the problem of 
finding stellar progenitors that have all of the requisite properties. 
SNe with evidence for circumstellar interaction represent ~ 10% of 
all core-collapse eve nts (,Smith et al.|2011[ l; ultraluminous SNe are 
much rarer still (e.g., |Quimby et al.|201 \) . This suggests that rather 
special stellar parameters are required to generate ~ Mq of ejecta 
in the last ~ year of stellar evolution. 

2 CONVECTION AND WAVE EXCITATION IN LATE 
STAGES OF STELLAR EVOLUTION 

Table [T] summarizes some of the key properties of core fu- 
sion and convection during the late stages of stellar evolution 
(based on |Woosley et al.|2002| and |Kippenhahn & Weigert|1990| l. 
For a stellar core with a mass ~ M©, a convective luminos- 
ity Lconv ~ Lfusion correspouds to a typical convective velocity 
V, ~ 10(L„„/10'Lo)'/3(p,/10'g cm-')-'/' km s^', where p, is 
the central density. The corresponding convective Mach numbers 
are given in Table [T| and are > 0.01 for Ne fusion and later. 
Although the parameters given in Table 1 are motivated by the 
core properties of evolved massive stars, similarly vigorous con- 
vection occurs during late stage shell burning. For example, |Amett| 
|& Mealanl ( |2011) find convective luminosities and Mach numbers 
comparable to those given in Table ^ for shell O fusion during the 
hour prior to core collapse. 

Figure [T| shows a mode propagation diagram for a 40 Mq, 
Z = lO""* metallicity model during core O fusion (at which point 



R ^ 1700 i^G, T'eff ^ 4000 K, and Lpho,o„ ^ 10^ Lq). The model was 
evolved using the MESA ID stellar evolution code ( |Paxton et al.| 
[20TT| l with no mass loss. Convective boundaries are determined by 
the Schwarzschild criterion, and hydrogen- and non-burning con- 
vection zones have overshooting of 1% of the local pressure scale 
height. This progenitor choice is somewhat arbitrary and thus un- 
likely to actually be the optimal progenitor for wave-driven mass 
loss. We include this model to provide a quantitative example of 
the propagation diagram and likely mode properties during very 
late stages of stellar evolution. For this particular stellar model, the 
core convection has a Mach number Alconv ~ 0.01 and locally car- 
ries a luminosity Lconv ~ 0.1 l^f^smn ~ iC^ Lq, with a comparable 
convective luminosity in the shell C fusion present at r ~ O.O3i?0. 

Vigorous convection transfers some fraction of the turbulent 
kinetic energy into waves, both sound waves and internal gravity 
waves. For the specific stellar model in Fig.[T| the core convective 
region is bordered by a convectively stable region into which the 
internal gravity waves excited by convection can propagate. This 
is important because the roughly incompressible convective mo- 
tions are much more efficient at directly exciting internal gravity 
waves than acoustic waves ( [Goldreich & Kumar|199"0| l. Quantita- 
tively, convection carrying a luminosity Lco„v with a Mach number 
Mcom excites an internal gravity wave luminosity of ( [Goldreich &| 
[Kumar] 1990t 

Lwave ^ >'*^conv Lconv ~ 1^ I lO^^'L. I \ 01 / 

This analytic estimate of the fraction of the convective energy trans- 
ferred to internal gravity waves is consistent (at the order of magni- 
tude level) with simulations of internal gravity wave excitation by 
solar convection ( [Rogers & Glatz maier 2005 2006). 

The frequencies and wavelengths of the internal gravity waves 
are somewhat more uncertain than the total wave luminosity given 
in equation [T] If the excitation is dominated by motions in the 
convection zone itself, the characteristic frequency associated with 
most of the wave power is the convective turnover frequency, 
(jj ~ 0)^ = T~' . Likewise, the horizontal spatial scale is set by the 
size of the convective eddies, so that the characteristic angular de- 
gree of the excited modes is f ~ where H is the size of the 
convective eddies and r is the radius where the excitation occurs 
(e.g., [Kumar et al.|1999| ). For the core (and 'thick' shell) convec- 
tion relevant here, this would correspond to modes having C ~ few. 

The uncertainty in these estimates is that much of the exci- 
tation may instead occur in a convective overshoot layer ( )Garcia[ 
[Lopez & Spruit|[T99T] [Rogers & Gla tzmaier 2005 , 2006) l. This 
increases both the frequency and angular degree t of the excited 
waves because the length-scale over which the excitation happens 
is the thickness of the overshoot layer. For evolved massive stel- 
lar models the thickness of the overshoot layer is particularly un- 
certain because the vigorous convection itself generates significant 
mixing and may substantially modify the structure of the star rela- 
tive to that predicted by ID models ( [Meakin & Amett[2()06{ [Amett] 
[& Meakin|[20I Ij l. It is likely that the convective overshoot layer 
is significantly thicker in this circumstance than in the sun (where 
it is < 0.05// at the interior convective-radiative transition; e.g., 
[Christensen-Dalsgaard et ar|[2011^ . This favors the excitation of 
lower ( waves with a) ~ uic- For our analytic estimates we will use 

<jj = a(jar=a Alconv TT ~ « Alconv S 1 (2) 
H 

where a ~ 1 - 10 is a dimensionless number that encapsulates the 
uncertainty associated with the wave excitation and S i is the ^ = 1 



© 0000 RAS, MNRAS 000, 000-000 



Mass Loss in the Last Year of Massive Stellar Evolution 3 




log r |ff e 



Figure 1. Propagation diagram for a 40 M© star with Z = 10 during core O fusion, showing the Brunt- Vaisala frequency N (blue, solid line) and the Lamb 
frequency St (for t = \ and 10; green, solid) on the left ordinate, plotted as a function of radius. The convective mach number Alconv (gray, dashed) is also 
shown on the right ordinate. The thin red lines show the convective turnover frequency in the core, and ten times that frequency, to demonstrate the range of 
internal gravity wave frequencies the core convection is likely to excite. The inner and outer radius of the tunneling cavity (r,„ and To,,,) for an f = 1 mode of 
frequency oj^ are marked on the top abcissa. 



Lamb frequency. The characteristic wave frequencies oj = w,. and 
(jj = 10 a>^ are shown with the thin red lines in Fig.[T] 

3 TUNNELING, TRAPPING, AND DAMPING 
3.1 Internal Gravity Waves 

Depending on their frequency and angular degree, the energy sup- 
plied to internal gravity waves can either remain trapped in the cen- 
tral part of the star or tunnel out to the stellar envelope in the form 
of propagating sound waves. We now estimate the conditions under 
which the latter occurs. 

For internal gravity waves that become evanescent at a radius 
and tunnel to a radius r„^f, the effective damping rate of the oth- 
erwise trapped mode due to leakage to large radii is given by jUnno| 
|etal.|1989| > 

Vg / Tin \ (Ll ri„ \ 

r,n \roi,i/ A<N) \r„„,/ 

where A' = ^(^ + 1) and (N) is the average Brunt-Vaisala fre- 
quency in the propagating region. The radii r,„ and are labeled 
in Fig.[T]for £ = I modes with a> = ai^. Physically, the timescale for 
energy to leak out implied by equation [5] (~ Tje'i;) is given by the 
group travel time across the propagating region (r,„/vj) divided by 
the tunneling probability, i.e., the fraction of the mode energy that 
tunnels through the barrier in a given group travel time. 

In addition to the tunneling captured by equation [3] the out- 
going gravity waves can also be partially reflected by the rapid 
changes in composition at shell boundaries (see the spikes in A'^ in 
Fig.[T](; this occurs if the wavelength of the waves is larger than the 
thickness of the region over which the composition changes. The 
latter is set by convective overshoot at the base of the shell convec- 
tion zones and is probably a few percent of a scale height. Given 
this, we estimate that lower frequency gravity waves with a) ~ oj^ 
are reasonably in the WKB limit in which equation |3] applies, but 
for higher frequency waves with u) ~ lOiWc, the tunneling may be 



suppressed by an additional factor of ~ 10 due to the compositional 
boundaries. 

Internal gravity waves of frequency oj (given by eq. [2]( and 
degree i excited at the outer edge of the core convection zone at 
radius ~ r begin to tunnel at a radius ~ 3r (based on the width 
of the convectively stable region at ~ 0.01/?o in Fig.[TJ. The radius 
where the tunneling ceases (rou,) is determined by where w > 5f, so 
that the waves become propagating sound waves. For massive stel- 
lar models at radii ~ 10"^ - 10/?o> we find that the Lamb frequency 
can be reasonably approximated as a power-law in radius 5f oc r"*, 
with b ~ 1.2 - 1.5. Using equation |2] for the characteristic wave 
frequency, this implies rout/Hn - 0.3 (cf Alconv) ''''A''* and thus 

/ i\2[A/b] + l 

rieak ~ 3-'^ (ttM,„„vA-') OJ. (4) 

The dominant damping mechanisms for the gravity wave en- 
ergy trapped in the core of the star are nonlinear damping and neu- 
trino damping. The neutrino damping rate is 

7,~10t.l,'~100t-J, (5) 

where is the thermal (Kelvin-Helmholz) time and the factor of 
10 in front of r^,' is due to the strong temperature dependence of the 
neutrino reactions (so that a given perturbation in temperature gives 
rise to a yet stronger perturbation in neutrino cooling). The last ex- 
pression in equation|5]is a consequence of the fact that the duration 
of nuclear burning (tnu^ in Table [TJ in the late stages of massive 
stellar evolution is only a factor of ~ 10 longer than the thermal 
time. For some of the g-modes of interest it is possible that driving 
due to fusion may be stronger than the neutrino damping (see, e.g., 
[Murphy et al.|20C)4^ . The interaction between convective excitation 
and such driving could be very interesting, but to be conservative, 
we do not explicitly include this driving in our estimates. 

During O and Ne fusion, equation |5] implies that y,, ~ 3 x 
10"* Hz ~ (jj^, for w,. ~ 3 X 10"' as shown in Fig.[l] Thus, 
neutrino damping is rather effective for modes of frequency ~ uj^- 
For comparison, taking M^om ~ 0.01, we find that yieak ^ 7v for 
modes with f = 1, if a > 3. For 1 = 2, the condition is more 
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prohibitive: hiigh frequency modes liaving a > 10 are required for 
leakage to dominate neutrino damping. This demonstrates that only 
the power in the lowest degree modes with £ ^ 1 - 2 is capable of 
efficiently leaking out into the stellar envelope. These low £ modes 
are, however, those that are expected to be excited by the large-scale 
convection in late stages of stellar evolution. 

The nonlinear damping of internal gravity waves is somewhat 
more difficult to quantify. A useful measure of the nonlinearity in 
the propagating region is the dimensionless parameter k,.^r (the ra- 
dial displacement relative to the radial wavelength): when k^^,. > 1, 
the waves can locally overturn the stratification leading to efficient 
nonlinear damping. Using conservation of energy flux, ky^,. for a 
traveling internal gravity wave of frequency cj can be rewritten as 



1/2 



(6) 



where F„avc = L„avc/47rr- is the local wave energy flux. Quantita- 
tively, using the stellar progenitor shown in Fig. [T] and assuming 
L„ave ~ 10*Lo, we find that k^^^ <; 1 for radii < 0.015 /?o for 
a) ~ ci)^, while k^^,. < 0.01 at all radii of interest for a) ~ IOlj^. 
Wave breaking is thus unlikely to significantly limit the energies 
attained by higher frequency internal gravity waves, which also tun- 
nel the most effectively (this corresponds to larger a in eq.|4](. How- 
ever, wave breaking may be important for waves with frequencies 
~ cl)^. Indeed, Meaki n & Arnett| ( |2006^ see some evidence for mix- 
ing induced by g-modes breaking in the convectively stable region 
separating O and C-burning shells in their simulations of late-stage 
burning in massive stars. 

3.2 Sound Waves 

The outward energy flux in sound waves created as described in 
^3.1| dissipates well before reaching the stellar surface. In partic- 
ular, we find that the dominant dissipation of the sound waves is 
via radiative diffusion. The sound waves damp when the radiative 
diffusion time ~ l/(k^x) (where is the radiative conductivity) is 
shorter than the group travel time of the modes across a given scale- 
height of the star ~ ///c.,. This condition can be rewritten as 



(7) 



where we have assumed that gas pressure is not much smaller than 
radiation pressure (which is reasonable for the conditions of in- 
terest). Fig. |2] shows Lrad (which is ~ Lndd even though the en- 
velope is convective) and Ljamp as a function of radius for modes 
with a; ~ 3 Wc in the 40 Mq, Z = 10"'' progenitor used in Fig. 
[T] The sound waves damp by radiative diffusion in the envelope at 
r~ 30- 100 

The dissipation of the outgoing sound waves will drive con- 
vection in the outer stellar envelope. However, the convection 
ceases to be efficient when the wave luminosity is larger than the 
maximum energy that subsonic convection can carry Lrnax.conv = 
An pel- So long as Lwave ^ Lfiddi this will always occur in- 
side the photosphere. For the specific model shown in Fig. |2] 
Lwave > Lconv.max outside ~ 500i?o- There is ~ 5Mq of mass ex- 
terior to this radius. The inability of convection or radiative diffu- 
sion to carry the outgoing wave energy implies that the wave power 
will inevitably drive a strong outflow. In particular, the mass out- 
flow rate induced is likely to be M ~ Anr-pc^ (evaluated where 
L„ave ~ Lconv.max. SO that convectlon cannot carry the energy; this 
corresponds to the sonic point of the outflow). For the stellar model 
used in Figures[T|and|2] we find M ~ IOMq yr"' and L„avc ^ Mv^^^. 
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Figure 2. Critical luminosities that determine the dissipation of the out- 
going wave power Lwavc and tlie abiiity of convection to carry the ther- 
malized energy. Ail quantities are for the 40Mq, Z = 10"'* stellar modei 
from Fig.^ Outgoing acoustic waves damp by radiative diffusion when the 
background radiative flux in the steifar model L,-ad exceeds the critical lu- 
minosity Ldamp (see eq. |7] w = 3oJc for Ldamp in this example). Subsonic 
convection cannot carry the thermalized wave power (L„ave = IO'L© here) 
when Lwave S Lconv.max = Anr^pc] . Taken together, these results imply that 
the super-Eddington wave power cannot be earned by either convection or 
radiative difl'usion in the stellar envelope at r > SOOSq, leading to strong 
mass toss. There is ~ 5Mq exterior to this radius in this exampie. 



(for Lwave ^ 10' Lo), so that the wave power is capable of driving a 
sustained outflow. 

4 DISCUSSION 

The total energy released during Ne and O fusion in the last year 
of massive stellar evolution is ~ 10" ergs. We have shown that 
a significant fraction of this energy ~ IQ'"*"*' ergs - i.e., Lwave ~ 
lO"*'"*" erg s"' - is converted into internal gravity waves via the 
Mach > 0.01 convection that accompanies the enormous nuclear 
and neutrino luminosities in the cores of massive stars. We have 
further argued that a large fraction of the power in low angular de- 
gree (£ ~ few) modes can tunnel through to the stellar envelope and 
become outgoing acoustic waves. 

If the low degree modes carry a significant fraction of the total 
wave power, the outgoing energy flux in acoustic waves will be sig- 
nificantly super-Eddington. The dissipation of the acoustic waves 
in the stellar envelope then drives convection that attempts to carry 
a super-Eddington power ~ lO'"'"'*' erg s"' ~ 10 - lOOLndd- Such 
convection becomes inefficient inside the stellar photosphere (Fig. 
[2| and thus the end result of this wave energy deposition is almost 
certainly substantial mass loss. Assuming that the outgoing wave 
power can be maintained for a reasonable fraction of the duration 
of Ne and O core fusion, up to 10'*'"''** ergs is deposited in the stellar 
envelope. If the unbound material moves at ~ 100 - 1000 km s"', 
the wave power can unbind ~ 10"^ - IOMq of material. 

Matter moving at ~ 300km s"' ejected in the year prior to 
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core-collapse will end up at ~ 100 AU when the star explodes. 
The mass-loss rates, ejecta mass, and the radial extent of the ejecta 
estimated here are consistent with the conditions required to ex- 
plain SNe powered by circumstellar interaction between the SN 
shock and surrounding stellar ejecta. In particular, at the low end 
of our estimated mass loss rates (~ 10"'' - 10"^ Mq yr"'), circum- 
stellar interaction can produce typical Type Iln supemovae ( |Kiewe| 
|et al.|20r2} . More extreme cases with ejecta masses approaching 
~ 0.1 - IMq in the last year are comparable to what is needed 
to power the most luminous SNe ever detected (e.g., SN 2006gy, 
[Smith et al.|2007| and the ultraluminous Ics, |Quimby et al.|2011^ 
as well as SNe powered by shock breakout in a dense circumstellar 
wind (e.g., |Ofek et al.|2010 t. Moreover, our results provide an ex- 
planation for what is otherwise a fine tuning problem: why should 
the star happen to lose a non-negligible fraction of its mass in the 
last year of stellar evolution (~ 10"*"^ of its lifetime!)? 

For red super-giant progenitors, the mechanism of mass loss 
proposed here operates primarily during core Ne and O fusion. Dur- 
ing core C fusion, the energy flux in internal gravity waves is typi- 
cally well below Eddington and so is unlikely to modify the stellar 
mass loss. During Si fusion, the wave luminosities are enormous, 
but the sound crossing time of a giant (~ few months) is much 
longer than the duration of the burning phase and so the star un- 
dergoes core-collapse before waves can reach the stellar surface. 
Nonetheless, the outgoing wave power created during Si fusion 
may substantially modify the structure of the progenitor star at radii 
< 30/?o (see Meaki n & Arnett|2006[ l. This could change how the 
supernova shock interacts with the surrounding star. In addition, 
in compact progenitors that have already lost their hydrogen en- 
velopes, wave excitation during Si fusion could help drive substan- 
tial mass loss in the day prior to core-collapse. Such stripped enve- 
lope progenitors are of particular interest in the context of under- 
standing whether circumstellar interaction powers hydrogen-poor 
ultra-luminous SNe (e.g., |Quimby et al.|20I \\ . 

In the well-studied case of SN 2006gy - whose lightcurve 
is consistent with ~ IOMq of ejecta in the 8 years prior to core- 
collapse l |Smith et al.|2010| > - there is also evidence for a compa- 
rable amount of ejecta ~ 10^ years earlier (via an IR light echo; 
[Miller et al.|2010p . The wave-driven mass loss mechanism we have 
proposed cannot work in its present form for this earlier mass loss 
episode. Likewise, it is unlikely to be relevant to typical luminous- 
blue variable outbursts. 

Our model for wave-driven mass loss requires that a reason- 
able fraction of the internal gravity wave power excited by convec- 
tion reside in low £ modes with frequencies somewhat larger than 
the characteristic convective turnover frequency. The reason for the 
former condition is that high £ gravity waves cannot tunnel through 
to the stellar surface; their energy is instead trapped in the interior, 
where it is ultimately thermalized and radiated via neutrinos. The 
latter condition is required since lower frequency gravity waves will 
likely break locally within the g-mode propagating region and in- 
duce mixing (see eq. |6j such mixing could be very interesting in 
its own right: e.g., by bringing fresh fuel down to higher temper- 
atures where it could in principle combust and power an eruption 
like those studied by Dessart et al. 120101). In their simulations of O 
shell fusion, |Meakin & Amett |2006 ) found significant wave power 
in ^ = 4 modes; since these simulations covered only a quadrant of 
the star, these were the lowest t modes they could simulate. This is 
consistent with our hypothesis that the convection in the late stages 
of stellar evolution will be particularly large-scale and thus will ex- 
cite low £ modes. 

It is unclear whether core convection or shell convection is 



likely to be the most important source of waves for driving mass 
loss. Core convection tends to be more energetic and is likely larger 
scale, favoring the excitation of the low I modes that tunnel most 
effectively. However, a countervailing consideration is that waves 
excited during shell convection have a smaller barrier to tunnel 
through. Further numerical work to calibrate the power-spectrum 
of wave excitation is clearly needed, particularly during core O and 
Ne fusion. 

The basic energetics of core convection and wave excitation 
that we have highlighted apply to all massive stars. It is, however, 
clear on observational grounds that not all massive stars undergo 
extreme mass loss in the year prior to core-collapse. Given that our 
mechanism requires that a significant fraction of the gravity wave 
power excited by convection is in relatively low £ modes and that 
the tunneling cavity not be too spatially extended, we suspect that 
the answer to this apparent tension is that only in certain progen- 
itors (with particular mass, metallicity, and/or rotation) are these 
conditions satisfied. However, determining exactly which progeni- 
tors these are will require additional work. For example, as empha- 
sized by [Amett & Meakin| ( |20I l) , the structure of the core of the 
star itself depends on the mixing induced by the internal gravity 
waves and so standard one-dimensional models may not be suffi- 
ciently accurate to address these questions. 

A second question that we have not addressed is how the star 
responds to the wave power deposited in the stellar envelope. It is 
possible that in some cases the stellar structure adjusts in such a 
way as to suppress the fraction of the wave power that can tunnel 
to large radii before much of the stellar envelope has been shed. 
This would limit the total amount of mass loss via wave energy 
deposition. This question will be explored in detail in future work. 
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